Possible Contexts of Use for In Silico Trials Methodologies: A Consensus-Based Review.

The term "In Silico Trial" indicates the use of computer modelling and simulation to evaluate the safety and efficacy of a medical product, whether a drug, a medical device, a diagnostic product or an advanced therapy medicinal product. Predictive models are positioned as new methodologies for the development and the regulatory evaluation of medical products. New methodologies are qualified by regulators such as FDA and EMA through formal processes, where a first step is the definition of the Context of Use (CoU), which is a concise description of how the new methodology is intended to be used in the development and regulatory assessment process. As In Silico Trials are a disruptively innovative class of new methodologies, it is important to have a list of possible CoUs highlighting potential applications for the development of the relative regulatory science. This review paper presents the result of a consensus process that took place in the InSilicoWorld Community of Practice, an online forum for experts in in silico medicine. The experts involved identified 46 descriptions of possible CoUs which were organised into a candidate taxonomy of nine CoU categories. Examples of 31 CoUs were identified in the available literature; the remaining 15 should, for now, be considered speculative.

Estimation of left ventricular stroke work based on a large cohort of healthy children.

Left ventricular stroke work is an important prognostic marker to analyze cardiac function. Standard values for children are, however, missing. For clinicians, standards can help to improve the treatment decision of heart failures. For engineers, they can help to optimize medical devices. In this study, we estimated the left ventricular stroke work for children based on modeled pressure-volume loops. A lumped parameter model was fitted to clinical data of 340 healthy children. Reference curves for standard values were created over age, weight, and height. Left ventricular volume was measured with 3D echocardiography, while maximal ventricular pressure was approximated with a regression model from the literature. For validation of this method, we used 18 measurements acquired by a conductance catheter in 11 patients. The method demonstrated a low absolute mean difference of 0.033 J (SD: 0.031 J) for stroke work between measurement and estimation, while the percentage error was 21.66 %. According to the resulting reference curves, left ventricular stroke work of newborns has a median of 0.06 J and increases to 1.15 J at the age of 18 years. Stroke work increases over weight and height in a similar trend. The percentile curves depict the distribution. We demonstrate how reference curves can be used for quantification of differences and comparison in patients.

Virtual implantations to transition from porcine to bovine animal models for a total artificial heart.

Realheart total artificial heart (TAH) is a novel, pulsatile, four-chamber total artificial heart which had been successfully tested acutely in a porcine animal model. However, the bovine model is better suited for long-term testing and thus an evaluation of how the design would fit the bovine anatomy was required. Virtual implantation is a method that enables a computer simulated implantation based on anatomical 3D-models created from computer tomography images. This method is used clinically, but not yet adopted for animal studies. Herein, we evaluated its suitability in the redesign of the outer dimensions and vessel connections of Realheart TAH to transition from the porcine to the bovine animal model. Virtual implantations in combination with bovine cadaver studies enabled a series of successful acute bovine implantations. Virtual implantations are a useful tool to replace the use of animals in early device development and refine subsequent necessary in vivo experiments. The next steps are to carry out human virtual implantations and cadaver studies to ensure the design is optimized for all stages of testing as well as the final recipient.

Crucial Aspects for Using Computational Fluid Dynamics as a Predictive Evaluation Tool for Blood Pumps.

The suitability of computational fluid dynamics (CFD) as a regulatory tool for safety assessment of medical devices is still limited: A lack of standardized validation and evaluation methods impairs the quantitative comparability and reliability of simulation studies, particularly regarding the assessment of hemocompatibility. This study investigated important aspects of validation and verification for three common turbulence modeling approaches (laminar, k-ω shear stress transport [SST] and stress-blended eddy simulation [SBES]) and three different mesh refinements. Simulation results for pressure head, characteristic velocity, and shear stress for the benchmark blood pump model of the Food and Drug Administration critical path initiative were compared with its published experimental results. For the highest mesh resolution, all three models predicted the hydraulic pump characteristics with a relative deviation averaged over six operating conditions below 6.1%. In addition, the SBES model showed an accurate agreement of the characteristic velocity field in the pump's diffusor region (relative error <2.9%), while the laminar and SST model calculated significantly elevated and deviating velocity amplitudes (>43.6%). The ability to quantify shear stress is fundamental for the prediction of blood damage. In this respect, this study demonstrated that: 1) a close agreement and validation of both pressure head and characteristic velocity was feasible and 2) the shear stress quantification demanded higher near-wall mesh resolutions, although such high resolutions were not required for the validation of only pressure heads or velocity. Hence, a mesh verification analysis for shear stresses may prove significant for the development of credible CFD blood damage predictions in the future.

Fluid-structure interaction of a pulsatile flow with an aortic valve model: A combined experimental and numerical study.

The complex fluid-structure interaction problem associated with the flow of blood through a heart valve with flexible leaflets is investigated both experimentally and numerically. In the experimental test rig, a pulse duplicator generates a pulsatile flow through a biomimetic rigid aortic root where a model of aortic valve with polymer flexible leaflets is implanted. High-speed recordings of the leaflets motion and particle image velocimetry measurements were performed together to investigate the valve kinematics and the dynamics of the flow. Large eddy simulations of the same configuration, based on a variant of the immersed boundary method, are also presented. A massively parallel unstructured finite-volume flow solver is coupled with a finite-element solid mechanics solver to predict the fluid-structure interaction between the unsteady flow and the valve. Detailed analysis of the dynamics of opening and closure of the valve are conducted, showing a good quantitative agreement between the experiment and the simulation regarding the global behavior, in spite of some differences regarding the individual dynamics of the valve leaflets. A multicycle analysis (over more than 20 cycles) enables to characterize the generation of turbulence downstream of the valve, showing similar flow features between the experiment and the simulation. The flow transitions to turbulence after peak systole, when the flow starts to decelerate. Fluctuations are observed in the wake of the valve, with maximum amplitude observed at the commissure side of the aorta. Overall, a very promising experiment-vs-simulation comparison is shown, demonstrating the potential of the numerical method.

Development of an In Vitro PIV Setup for Preliminary Investigation of the Effects of Aortic Compliance on Flow Patterns and Hemodynamics.

The aorta with its compliance plays a major role in hemodynamics as it saves a portion of ejected blood during systole which is then released in diastole. The aortic compliance decreases with increasing age, which is related to several cardiovascular imparities and diseases. Changes in flow patterns and pressure curves, due to varying aortic compliance, are difficult to investigate in vivo. As a result, the aim of the present work was to develop an in vitro setup enabling standardized investigations on the effect of compliance changes on flow patterns and pressure curves. Therefore an experimental setup with an anatomically correct silicone phantom of the aortic arch was developed, suitable for optical flow measurements under pulsatile inflow conditions. The setup was developed for precise adjustments of different compliances and optical flow measurements. Particle image velocimetry measurements were carried out downstream of the aortic valve in the center plane perpendicular to the valve with compliance adjusted between 0.62 × 10-3 to 1.82 × 10-3 mmHg-1. Preliminary results of the in vitro investigations showed that decreases in compliance results in significant increases in pressure changes with respect to time (dp/dt) and altered pressure curves in the aortic arch. In terms of flow, an increased aortic stiffness lead to higher mean velocities and decreased vortex development in the aortic sinuses. As in vivo validation and translation remains difficult, the results have to be considered as preliminary in vitro insights into the mechanisms of (age-related) compliance changes.

Fluid Dynamics in Rotary Piston Blood Pumps.

Mechanical circulatory support can maintain a sufficient blood circulation if the native heart is failing. The first implantable devices were displacement pumps with membranes. They were able to provide a sufficient blood flow, yet, were limited because of size and low durability. Rotary pumps have resolved these technical drawbacks, enabled a growing number of mechanical circulatory support therapy and a safer application. However, clinical complications like gastrointestinal bleeding, aortic insufficiency, thromboembolic complications, and impaired renal function are observed with their application. This is traced back to their working principle with attenuated or non-pulsatile flow and high shear stress. Rotary piston pumps potentially merge the benefits of available pump types and seem to avoid their complications. However, a profound assessment and their development requires the knowledge of the flow characteristics. This study aimed at their investigation. A functional model was manufactured and investigated with particle image velocimetry. Furthermore, a fluid-structure interaction computational simulation was established to extend the laboratory capabilities. The numerical results precisely converged with the laboratory measurements. Thus, the in silico model enabled the investigation of relevant areas like gap flows that were hardly feasible with laboratory means. Moreover, an economic method for the investigation of design variations was established.

Using 3D Physical Modeling to Plan Surgical Corrections of Complex Congenital Heart Defects.

Background Understanding the anatomy and physiology of congenital heart defects is crucial for planning interventions in these patients. Congenital heart procedures often involve complex three-dimensional (3D) reconstructions. Excellent imaging techniques are required to depict all anatomical details. We have used and evaluated fast 3D prototyping technology for reconstruction and planning of corrections of complex congenital heart defects. Materials and Methods 3D physical models were constructed from contrast-enhanced computed tomography (CT) datasets of patients with complex congenital heart defect. Two different commercially available printing technologies were used and their clinical application compared. Results Physical models of three different patients were used for preoperative surgical planning. All models showed good correspondence to patient anatomy. Both printing technologies gave excellent results. Conclusion Physical models could be easily constructed with the use of CT datasets. The printing process could be done efficiently, quite rapidly, and cost effectively. Surgical corrections could be planned based on these models.

Hemodynamic analysis of outflow grafting positions of a ventricular assist device using closed-loop multiscale CFD simulations: Preliminary results.

Subclavian arteries are a possible alternate location for left ventricular assist device (LVAD) outflow grafts due to easier surgical access and application in high risk patients. As vascular blood flow mechanics strongly influence the clinical outcome, insights into the hemodynamics during LVAD support can be used to evaluate different grafting locations. In this study, the feasibility of left and right subclavian artery (SA) grafting was investigated for the HeartWare HVAD with a numerical multiscale model. A 3-D CFD model of the aortic arch was coupled to a lumped parameter model of the cardiovascular system under LVAD support. Grafts in the left and right SA were placed at three different anastomoses angles (90°, 60° and 30°). Additionally, standard grafting of the ascending and descending aorta was modelled. Full support LVAD (5l/min) and partial support LVAD (3l/min) in co-pulsation and counter-pulsation mode were analysed. The grafting positions were investigated regarding coronary and cerebral perfusion. Furthermore, the influence of the anastomosis angle on wall shear stress (WSS) was evaluated. Grafting of left or right subclavian arteries has similar hemodynamic performance in comparison to standard cannula positions. Angularity change of the graft anastomosis from 90° to 30° slightly increases the coronary and cerebral blood flow by 6-9% while significantly reduces the WSS by 35%. Cannulation of the SA is a feasible anastomosis location for the HVAD in the investigated vessel geometry.

Investigation of hemodynamics during cardiopulmonary bypass: A multiscale multiphysics fluid-structure-interaction study.

Neurological complications often occur during cardiopulmonary bypass (CPB). Hypoperfusion of brain tissue due to diminished cerebral autoregulation (CA) and thromboembolism from atherosclerotic plaque reduce the cerebral oxygen supply and increase the risk of perioperative stroke. To improve the outcome of cardiac surgeries, patient-specific computational fluid dynamic (CFD) models can be used to investigate the blood flow during CPB. In this study, we establish a computational model of CPB which includes cerebral autoregulation and movement of aortic walls on the basis of in vivo measurements. First, the Baroreflex mechanism, which plays a leading role in CA, is represented with a 0-D control circuit and coupled to the 3-D domain with differential equations as boundary conditions. Additionally a two-way coupled fluid-structure interaction (FSI) model with CA is set up. The wall shear stress (WSS) distribution is computed for the whole FSI domain and a comparison to rigid wall CFD is made. Constant flow and pulsatile flow CPB is considered. Rigid wall CFD delivers higher wall shear stress values than FSI simulations, especially during pulsatile perfusion. The flow rates through the supraaortic vessels are almost not affected, if considered as percentages of total cannula output. The developed multiphysic multiscale framework allows deeper insights into the underlying mechanisms during CPB on a patient-specific basis.

A numerical framework to investigate hemodynamics during endovascular mechanical recanalization in acute stroke.

Ischemic stroke, caused by embolism of cerebral vessels, inflicts high morbidity and mortality. Endovascular aspiration of the blood clot is an interventional technique for the recanalization of the occluded arteries. However, the hemodynamics in the Circle of Willis (CoW) are not completely understood, which results in medical misjudgment and complications during surgeries. In this study we establish a multiscale description of cerebral hemodynamics during aspiration thrombectomy. First, the CoW is modeled as a 1D pipe network on the basis of computed tomography angiography (CTA) scans. Afterwards, a vascular occlusion is placed in the middle cerebral artery and the relevant section of the CoW is transferred to a 3D computational fluid dynamic (CFD) domain. A suction catheter in different positions is included in the CFD simulations. The boundary conditions of the 3D domain are taken from the 1D domain to ensure system coupling. A Eulerian-Eulerian multiphase simulation describes the process of thrombus aspiration. The physiological blood flow in the 1D and 3D domains is validated with literature data. Further on, it is proved that domain reduction and pressure coupling at the boundaries are an appropriate method to reduce computational costs. Future work will apply the developed framework to various clinical questions.

Effect of pulmonary conduit oversizing on hemodynamics in children.

PURPOSE: Implanting the largest valved conduit possible - oversizing - to reconstruct an absent connection from the right ventricle to the pulmonary artery in certain types of congenital heart defects has been suggested as a compensating measure for somatic outgrowth of the patient. However, one effect that has not been investigated yet is the hemodynamic consequence. For this purpose, virtual implantation and flow simulations were conducted in this study. METHODS: To isolate the effects of conduit oversizing on the hemodynamics observed after conduit implantation and outgrowth, calculated wall shear stresses (WSS) of image-based computational fluid dynamic (CFD) simulations were used as indicators. Three different sizes of valved conduits (20 mm, 22 mm, and 24 mm), including the largest possible conduit size, virtually implanted in a child-sized healthy pulmonary artery and the corresponding adult-sized model were investigated. RESULTS: The child and adult models show a decrease of the mean WSS (approx. 26%) in the whole domain with an increase of the conduit size. When looking at the mean WSS at the anastomosis, for the child model the WSS is significantly increased (approx. 40%) when oversizing (Z-score +3.21). In contrast, the stresses are decreased for the adult model (34%) when using the largest conduit (Z-score +0.25). CONCLUSIONS: Based on the results of this study, it must be considered that choosing a prosthesis size that will lead to high WSS and an associated intimal reaction, possibly leading to stenosis, can defeat the benefit of having a nominally larger orifice area directly after implantation.

In vitro flow investigations in the aortic arch during cardiopulmonary bypass with stereo-PIV.

The cardiopulmonary bypass is related to complications like stroke or hypoxia. The cannula jet is suspected to be one reason for these complications, due to the sandblast effect on the vessel wall. Several in silico and in vitro studies investigated the underlying mechanisms, but the applied experimental flow measurement techniques were not able to address the highly three-dimensional flow character with a satisfying resolution. In this work in vitro flow measurements in a cannulated and a non-cannulated aortic silicone model are presented. Stereo particle image velocimetry measurements in multiple planes were carried out. By assembling the data of the different measurement planes, quasi 3D velocity fields with a resolution of~1.5×1.5×2.5 mm(3) were obtained. The resulting velocity fields have been compared regarding magnitude, streamlines and vorticity. The presented method shows to be a suitable in vitro technique to measure and address the three-dimensional aortic CPB cannula flow with a high temporal and spatial resolution.

Numerical washout study of a pulsatile total artificial heart.

PURPOSE: For blood pumps with long term indication, blood stagnation can result in excessive thromboembolic risks for patients. This study numerically investigates the washout performance of the left pump chamber of a pulsatile total artificial heart (TAH) as well as the sensitivity of the rotational orientation of the inlet bileaflet mechanical heart valve (MHV) on blood stagnation. METHODS: To quantitatively evaluate the washout efficiency, a fluid-structure interaction (FSI) simulation of the artificial heart pumping process was combined with a blood washout model. Four geometries with different orientations (0°, 45°, 90° and 135°) of the inlet valve were compared with respect to washout performance. RESULTS: The calculated flow field showed a high level of agreement with particle image velocimetry (PIV) measurements. Almost complete washout was achievable after three ejection phases. Remains of old blood in relation to the chamber volume was below 0.6% for all configurations and were mainly detected opposite to the inlet and outlet port at the square edge where the membrane and the pump chamber are connected. Only a small variation in the washout efficiency and the general flow field was observed. An orientation of 0° showed minor advantages with respect to blood stagnation and recirculation. CONCLUSIONS: Bileaflet MHVs were demonstrated to be only slightly sensitive to rotation regarding the washout performance of the TAH. The proposed numerical washout model proved to be an adequate tool to quantitatively compare different configurations and designs of the artificial organ regarding the potential for blood stagnation where experimental measurements are limited.

Combined computational and experimental approach to improve the assessment of mitral regurgitation by echocardiography.

Mitral regurgitation (MR) is one of the most frequent valvular heart diseases. To assess MR severity, color Doppler imaging (CDI) is the clinical standard. However, inadequate reliability, poor reproducibility and heavy user-dependence are known limitations. A novel approach combining computational and experimental methods is currently under development aiming to improve the quantification. A flow chamber for a circulatory flow loop was developed. Three different orifices were used to mimic variations of MR. The flow field was recorded simultaneously by a 2D Doppler ultrasound transducer and Particle Image Velocimetry (PIV). Computational Fluid Dynamics (CFD) simulations were conducted using the same geometry and boundary conditions. The resulting computed velocity field was used to simulate synthetic Doppler signals. Comparison between PIV and CFD shows a high level of agreement. The simulated CDI exhibits the same characteristics as the recorded color Doppler images. The feasibility of the proposed combination of experimental and computational methods for the investigation of MR is shown and the numerical methods are successfully validated against the experiments. Furthermore, it is discussed how the approach can be used in the long run as a platform to improve the assessment of MR quantification.

Implementation of intrinsic lumped parameter modeling into computational fluid dynamics studies of cardiopulmonary bypass.

Stroke and cerebral hypoxia are among the main complications during cardiopulmonary bypass (CPB). The two main reasons for these complications are the cannula jet, due to altered flow conditions and the sandblast effect, and impaired cerebral autoregulation which often occurs in the elderly. The effect of autoregulation has so far mainly been modeled using lumped parameter modeling, while Computational Fluid Dynamics (CFD) has been applied to analyze flow conditions during CPB. In this study, we combine both modeling techniques to analyze the effect of lumped parameter modeling on blood flow during CPB. Additionally, cerebral autoregulation is implemented using the Baroreflex, which adapts the cerebrovascular resistance and compliance based on the cerebral perfusion pressure. The results show that while a combination of CFD and lumped parameter modeling without autoregulation delivers feasible results for physiological flow conditions, it overestimates the loss of cerebral blood flow during CPB. This is counteracted by the Baroreflex, which restores the cerebral blood flow to native levels. However, the cerebral blood flow during CPB is typically reduced by 10-20% in the clinic. This indicates that either the Baroreflex is not fully functional during CPB, or that the target value for the Baroreflex is not a full native cerebral blood flow, but the plateau phase of cerebral autoregulation, which starts at approximately 80% of native flow.